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In recent years a number of white dwarfs has been observed with very high surface 



' ,'. would be much higher ( ~ 10^^ G). In this paper, we analytically study the effect of 

^0 \ high magnetic field on relativistic cold electron, and hence its effect on the stability 

^ • and the mass-radius relation of a magnetic white dwarf. In strong magnetic fields, 

O , the equation of state of the Fermi gas is modified and Landau quantization comes 

^ I into play. For relatively very high magnetic fields (with respect to the energy density 
of matter) the number of Landau levels is restricted to one or two. We analyse the 

■ equation of states for magnetized electron degenerate gas analytically and attempt 

m ■ 

^2 . to understand the conditions in which transitions from the zero-th Landau level to 
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first Landau level occur. We also find the effect of the strong magnetic field on the 

■ star collapsing to a white dwarf, and the mass-radius relation of the resulting star. 
We obtain an interesting theoretical result that it is possible to have white dwarfs 
with mass more than the mass set by Chandrasekhar limit. 
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I. INTRODUCTION 



The origin of high magnetic fields in compact stars is explained by the fossil field hypoth- 
esis, as proposed by Ginzburg [ij and Woltjer 2]. The magnetic fiux 0b ~ AnBR^ of a star 
is conserved during its evolution, thus a degenerate collapsed star is expected to have a very 
high magnetic field if the original star had a magnetic field, say, B ~ 10^ G 3| j4|. Neutron 
stars were of the main interest in this category, but in recent years some white dwarfs have 
been found to have quite high surface magnetic fields with field strength varying from 10^ G 
to 10^ G js-S]. It is highly intuitive that the fields near the core of the white dwarf would 
be much higher. Ostriker and Hartwick jo] constructed models of white dwarf with central 
magnetic field 10^^ G but a much smaller field in surface. The maximum limit of the field 
strength is set by the virial theorem 



2T + W + 3U + M = 0, 



(1) 



where T is the total kinetic energy, W the gravitational potential energy, U the total internal 
energy and Ai the magnetic energy. Since T and U are both positive definite, the maximum 
magnetic energy is always less than the total gravitational energy in equilibrium. For a star 
of mass M and radius R this gives 



Bqnax ~ 



2 X 10* 



M 



(2) 



where Mq and Rq respectively denote solar mass and radius. For a white dwarf this sets 
the limit of i? ~ 10^^ G at the center but relatively lower fields outside. 

The mass-radius relation for a non-magnetic relativistic white dwarf as determined by 
Chandrasekhar [icj sets its maximum mass to 1.44 Mq such that the electron degeneracy 
pressure is adequate for counterbalancing the gravitational collapse. This mass limit is 
strengthened by higher central densities j3|. The effect of weak magnetic field B < 10^^ G 
on white dwarfs was studied by Suh et al. [ll| by applying Euler-MacLaurin expansion on 
the equation of state for a degenerate electron gas in a magnetic field 12[. They found that 
both the mass and radius of the white dwarf increase in presence of magnetic fields. 

In this report we analyse the case of relativistic white dwarfs with magnetic fields 
B > 10^^ G. In order to have a white dwarf of such magnetic field, the original solar 
type star should have a magnetic field ~ 10^ G, from the fiux freezing theorem. Existence 
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of such stars is not ruled out j^. First we analyse the effect of such strong magnetic field 
on the equation of state of degenerate electron gas at zero temperature and analyse some 
of its properties, subsequently we consider considerable higher magnetic field for the same 
and then study explicitly the mass-radius relation for the case of just one Landau level 
(produced by high magnetic field). The results of higher Landau level occupancy (hence 
lower magnetic field) are addressed quantitatively from the equation of state. In doing so, 
we find for central densities p ~ 10^" fi'/cc the maximum mass of the white dwarf can be 
much greater than the Chandrasekhar limit M ~ IAAMq. 



II. EQUATION OF STATE IN HIGH MAGNETIC FIELD 



A. Density of states 

The energy eigenstates of the free electrons in the magnetic field are quantized to what 
is known as Landau orbitals. The electrons motion perpendicular to the magnetic field is 
no longer independent but is quantized. This can be seen by solving the time-independent 
Schrodinger equation of a particle in a magnetic field B directed along z— axis given by [l3| 



{p, + eBy/cf ^ pI ^ pI 



^ - —B-^ = E^, (3) 



2me 2me 2m 

where Px,Py,Pz are the components of linear momentum, y is the arbitrary position on the 
y— axis, me and e are the mass and charge of the electron respectively, c the speed of light, 
fi the magnetic moment, s the magnitude of spin, a = E the energy eigenvalue. We 
obtain the classical quantization of energy levels of the system given by 

= + Sr,- 

where uh = -^^^ is the critical cyclotron frequencies at which quantization occurs, h the 
Planck's constant, and u = (/ + 1 + a), gives the Landau levels where / is the principal quan- 
tum number for the electron. Relativistically we obtain energy by solving Dirac equation in 
a magnetic field, hence the energy above is modified to 

E^ = ^Jmlc^ + pIc^ + 2venBc. (5) 

Thus electrons with up (f) spin and down (4) spin have different energies. We see that the 
ground state [y = 0) has degeneracy 1, while the Landau levels from z/ = 1 have degeneracy 



4 



2. The spin J, electrons can occupy Landau levels v = 0,1,2,3..., while the spin ones 



take Landau levels u = 1,2,3... IJ]. Note that the motion of the electrons in the x — y 
plane is coupled. Therefore the density of states of the electron will change, as the motion 
is restricted and quantized in the plane. As the motion in the x — y plane is quantized, the 
phase space occupied by the electrons will change. The number of states per unit volume in 
the interval Ap^ for a given Landau level u is (^fj^) Apz- The modified density of state is 
then 

2 p ^ 2—}_^g,odp.iu), (6) 

where g^o = (2 — ^o.j/)? is the degeneracy in each Landau level, and dpzi^u) is the small 
element of component of momentum in the 2— direction in the i/— th Landau level. Thus 
each Landau level has its own distribution of states. The separation between two Landau 
levels depends on the strength of the magnetic field. For high magnetic fields the separation 
of the Landau levels is large, hence electrons with low energy (non-relativistic) can only 
occupy the ground state. As the strength of the magnetic field decreases, the separation 
between the levels decreases. Hence it becomes energetically favourable for the electrons to 
jump to higher levels, thus the number of occupied Landau level increases. Similarly in the 
case of relativistic electrons, if the magnetic field is low, the separation of the Landau levels 
is comparable to the rest mass energy of the electrons, and hence the electrons can freely 
move between the Landau levels, thus it will behave as continuum. But again in presence of 
very high magnetic fields, the electrons, although relativistic, cannot jump to higher levels. 
The electrons can become relativistic in two different cases: (1) when the density is high 
enough such that the mean Fermi energy of the electron exceeds rest mass energy of the 
electron, (2) when the term in the energy associated with the cyclotron frequency of the 
electron exceeds the rest mass of the electron. 



B. Equation of state 

We now analyse the equation of state for the relativistic electrons in magnetic field, as 
12| . We define the Fermi momentum, analogous to that in the non- 



was obtained earlier 



magnetic case, as p]? = pi + ^i^^^. Therefore for u-th Landau level the 2;— component of 



the momentum is given by Pz{i^) = \l Pf ~ ^'^'^^^ ■ We introduce a convenient dimensionless 



2 3 

parameter xp = as relativity parameter, and = ^^^^^^ = 4.14 x 10^'^ G, as the critical 
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magnetic field giving rise to the significant effect due to the Landau quantization. We further 
define 7 = -f-, and A = the Compton wavelength of the electron. Using these definitions 
we can write Pz(z/) as 

x^{u) = x{u) = - 2i/7, (7) 



Pz 



where x{v) denotes the relativity parameter of the 2;— component of momentum for z/-th 
Landau level. The number density is then given by 

where I'max is the maximum number of Landau level that can be filled which is determined 
by the condition that ^^(z/) is real. Hence 



< i^max = Integer ( 



,Xf 



Integer 



.62-1, 



27' ° ' 27 ' 

where ep = -[^^ which is the dimensionless chemical potential and from equation ([5]) 
1 + Xp. This is related to matter density by 



(9) 



P = PemnUe = n^Q, 



(10) 



where He is the mean molecular weight given by A/Z, Z the atomic number, A the mass 
number, and m„ the mass of the neutron. Note that = yUg^n which has the dimension of 
mass and denoting the effective mass per electron. Its value depends on the constituents of 
the white dwarf. In our case value of /le is taken to be 2, and hence B = 2m„. The electron 
energy density at zero temperature is given by 12| 

Xp 

E = > ■ Q,M I E,Ax(u), (11) 



(27r)2A3 



•^max p 

giyo / Ei,dx{u) 

u=0 { 



and hence 



u=0 



(27r)2A3 



,2 r 



x(z/)ei. + (1 + 2Z/7) ln[ 



x{iy) + e^, 
VI + 2z/7 



(12) 



Then the pressure of the Fermi gas at zero temperature can be found from the relation 

P = J2^Fx{iy) - E, (13) 

which gives the relation 

rxM + ep- 



'jguomeC^ 



u=Q 



(27r)2A3 



x{iy)eF - (1 + 2Z/7) ln[ 



(14) 
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III. ONE LEVEL AND THE TWO LEVEL SYSTEMS AND EQUATION OF 

STATES 

First we consider the one level system in which only the ground level is occupied. This 
is possible at a very high magnetic field. Expanding equation (jH]) to first term (for z/ = 0) 
and with equation f fTOj) we obtain 

from which we further obtain 

XF = P = (16) 

7/iem„ K 

where K = Pef^n (2^)2 • Here we obtain Xi;' in terms of density and hence the relation of 
Fermi energy to the density of state. As we increase the density of the matter, we in turn 
increase the Fermi energy associated with it, hence electrons acquire probability to jump 
to higher Landau levels and therefore v can not be set 0. From equation (j9]) we find that 
once the Fermi energy and v are fixed, the magnetic field and the maximum density are 
restricted. The equation of state f|T^ then reduces to 



for u = 0. 

For the two level system (when z/ = 0, 1) 



u=0 

and the pressure then reduces to 



= ^ E ((2 - ^0,.) \/4 - 27^) = ^ ( + 2\/4 - 27 ) (18) 



P = ' , xpep + 2eF\ xl - 27 - ln(x^ + e^) - 2(1 + 27) ln( 



(27r)2A3 \^ - - ■ - V ^ ' V - ■ y ' u y .JYT^ 

The density can be written in terms of Xi;' as 



(19) 



p = KyxF^2^x\-2^y (20) 

Hence, 

-p + 2v/6ir27 + p2 
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which is always positive definite. The other solution is neglected as infeasible solution as 
density cannot be negative for real particles. 

Substituting this Xi? in equation (|T9i) we obtain 



^ (tp (sKx, -p + + P=) - e(K + 21(1) In [ "V^, 



where 

Xi 

and 



_ 2 

n2 



'--t^^ — — '- (^^' 

We must be careful that the two level system is only valid above a density given by equation 
f l28|) . The net equation of state is represented by two functions given by: 

Equation (16) for p < pt (z^ = 0) (25) 
Equation (21) for p > Pt (^^ = !)• 

where pt is the density at which transition occurs between the Landau levels. A diagrammatic 
representation of the equations ( |T71) and ( |T9l) is given in FiglH 

Here the values of K will depend on the limiting magnetic field taken. As for one level 
system the limiting magnetic field is much higher than that of two level system, the values of 
K will be in general different between one level and two level systems. Hence, the steepness 
of the equation of state representing ground Landau level will be different between the 
systems restricted to one level and that restricted to higher levels. The higher level systems 
can be computed in a similar algebraic way. It must be noted that the value of 7 obtained 
from equation ([9]) corresponds to the lower limit for the given value oi ep- For example, 
for say maximum Fermi energy Epmax = '^OrrieC^, if 7 = 199.5 it corresponds to one level 
system, but for any other values greater than 199.5 will not alter the condition for the same 
level system. 
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FIG. 1: Equations of state for one level (solid line) and two level (dashed line) systems as rep- 
resented by equations p!7j) and ([19]) for ej? = 20. The kink in dashed curve at pt is found from 
equation (pO|) . 



The electrons in the system become relativistic when the factor ep becomes more than 
one. For relatively low magnetic fields, the spacing between Landau levels becomes small 
and Vmax becomes large. Thus relativistic electrons can freely jump between the states as 
said earlier. In this case the levels can be treated as continuum and the summation over 
states can be replaced by an integral and the resulting equation gives back the normal non- 
magnetic relativistic equation of state P = K'p^ with F = 4/3, just the case which was 
analysed by Chandrasekhar [l^. Note that F = 4/3 corresponds to the equation of state 
of relativistic but non-magnetized electrons giving rise to the Chandrasekhar mass limit, 
which we recover in a particular range of density even in the magnetized system. In just the 
opposite case, which we are interested in, we consider just one/two Landau level(s), i.e. the 
ground/first excited state, is filled in. From equation ([9]) we fix Umax to a value less than 1 
(2) such that only the ground level (and first level) is (are) occupied. The electrons in the 
ground level would all be in the spin J, state. From equation ([9]) we can see that once we fix 



A. Discussion of equation of states 
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i^max and Xp, then the minimum magnetic field is automatically fixed. 

If there are finite Landau levels, then the ground level would be filled first. If the density 
is increased further, then the ground state will saturate (it cannot contain more electrons) 
and there will be a phase transition like phenomena, in which electrons would start filling 
the first excited Landau level. This process will continue to higher Landau level transitions. 
During these phase transition like phenomena, the E—p relation shows a sharp discontinuity. 
Thus ^ = oo. An analytic condition can be found out for the critical densities at which 



this transition occurs by differentiating equation (fT2|) with respect to p. In for two states, 
we note there is only one unique denominator which is not positive definite, hence can be 
set to zero which makes ^ = 00. Thus this gives us the condition for the critical density. 
Solving the equation 

27 = xl (26) 
and equation (|2T|) simultaneously we obtain 

pt = y^K. (27) 

This gives the critical density below which the electrons cannot occupy the first Landau level 
(z/ = 1) and would all be in the ground level (z/ = 0) and hence it would be inappropriate 
to use the two level equation of state. This result can also be estimated by using ^/2r/ < xp 
which gives on combining with xp = -j^ the same formula as equation 0271) . Similarly, the 
transition from first level to second level would then be given by the equation 

-p + 2^QK^^ + p^ 



V2 X 27 < ^ ' " ■ (28) 

The critical density for the transition from the first to second levels is given by 

Pt = 2(1 + V2)^K. (29) 

The higher phase change also occurs under the same principle. Therefore we can write the 
formula of critical density for the transition from [u — l)-th level to z/-th level for 1 < u < 2 

as 

p,(z/) = (v^ + 2v/2(z/-l))V^ir, (30) 
where it must be remembered that K is a. function of 7. 
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IV. DETERMINING MASS AND RADIUS OF THE WHITE DWARF 

ANALYTICALLY 

For a magnetized star in hydrostatic equilibrium, we require to solve the condition for 
equilibrium given by Q 

where P is the electron degeneracy pressure, varying with the radial coordinate of the star 
r, playing the role to balance the gravitational force, G the Newton's gravitation constant. 
Now we choose the equation of state as a polytropic relation 

P = K'p^, (32) 

where K' is a dimensional constant, 

P = Pee"", (33) 

and 

r = at (34) 
where Pc is the central density of the white dwarf and a is a constant defined by 



(n + 1)K' 
AnG 



2 (1-n) 

Pc'" . (35) 



With the use of equations (1321) . (1331) . (l34l) and (!35ll . equation (l3Til reduces to famous Lane- 
Emden equation 

where F = 1 + The Lane-Emden equation can be solved for a given n using boundary 
conditions 

e{i = 0) = 1 (37) 

and 

For n < 5, the value of 9 falls to zero at a finite ^ say which basically defines the surface 
of the star where pressure goes to zero. The physical radius is given by 

R = a^o. (39) 
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Here we notice that n > — 1 so that a (see equation (135|) ) must be real and so does R. 

Since the equation of state in presence of magnetic field cannot be written in a simple 
polytropic form as in equation ( l32l) . we fit equation IHM for various density ranges. The 
actual equation of state is reconstructed by using multiple fits in various density range. The 
values for K' and F in various ranges of density, which are found from the fitting function, 
also carry information about the magnetic field of the system. The idea behind the fit is 
to be able to solve the Lame-Emden equation piecewise ranges of density and to obtain an 
idea of the mass-radius relation from it easily. The mass of the star is then given by 

M = j Airpr'^dr = Airp^a^ ^'^O^'d^ (40) 

Now we consider only one level system explicitly, hence fit the one level equation of state 
and obtain different values of K' and F. We see from Fig. [1] that the equation of state 
representing the one level system for e^? = 20 is given by 

P = 2 X 10^ (|pV34 X 1016 + p2 _ 34 X 10^*^ In (^{p + ^34 x lO^^ + p2) iq x iq-w'^ ^ , (41) 

which is obtained from equation (fT7|) . This curve is fitted with two functions 

P(p) = 6.98 X lO-V^'*^^ for p < 1.4 X 10^ 
P(p) = 1.98 X lOV for P > 1-4 X 10^ 

shown by Fig. EJ^b). Here we consider a uniform magnetic field. In very low densities where 
we consider the pressure goes to zero, the equation of state defined above cannot be used as 
it becomes non-linear in nature. We have extended the fitting curve until this region arises. 
However, this region is a very small part of the total density region, hence would not affect 
our results effectively. Similarly, the fitting of the one level equation of state for ei;' = 2 is 
given by Fig. |2]^a). Then we integrate equation ( HOl) in two ranges of density to obtain 

M = Anp, (^aimf ^O^d^ + «(^2)'(^^° ^'O^'d^ - ^'^"^^o) , (42) 

which gives 

M = 47rp, [a{n,)%'e'{^,)\ + ain^fH^o'e'm - \^^'e'{m) , (43) 
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FIG. 2: Equations of state for one Landau level system, when (a) ep = 2, (b) ep = 20. Solid and 
broken lines respectively indicate original and fitting results, both overlap each other perfectly. 
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where is the radius corresponding to the intermediate density where the T changes its value 
drastically in the equation of state. We then substitute this intermediate density (p(^i)) 
in equation ( l33l) with a particular value of central density and corresponding polytropic 
exponent and find the corresponding 6. From the solution to Lane-Emden equation we find 
the corresponding values of ^ and .^^|^'(^)|. The radius in this case will be 

R = aini)^i + ain2)i^o-^i), (44) 

where the values of a(ni) and 0(^2) are obtained from equation (!35i) by substituting the 
corresponding values of n and K' . 

V. MASS-RADIUS RELATION 

We now discuss the mass-radius relation for the highly magnetized star. We again consider 
to study the extreme case having only one Landau level. For the higher level occupancy, as 
shown by the equation of state in Fig. [1] by the dashed curve for a two level system, the 
magnetic effect will be softened and we will discuss the results qualitatively. From equation 
fH5]) we obtain mass as a function of central density and equation flH|) gives the radius of 
the star as a function of central density, for one level systems. The values of coefficients 
are calculated from the values of K' and F from the fitting curves in various zones and 7 
is determined from the Fermi energy. We can eliminate the central density from equations 
(143|) and (jH]) and obtain mass as a function of radius. In order to express the results in 
convenient units we write the mass in the units of solar mass and radius in the units of 
10^ cm. The values of K' and F obtained from the fitting are used in equations fl43p and 
(H4|) which gives for ep =20 

M/Mq = 3.81291 X 10""pc + 2.86571 x lO^^^c'^^ (45) 

and 

= 0.05839 + 3.196 x 10"Vc'^®- (46) 

Now eliminating pc from above equations we obtain the mass-radius relation. Figure [3] shows 
the mass-radius relation for ep = 2,20. The extreme right point of the curves (maximum 
radius) denotes the typical maximum central density of white dwarfs for the respective values 
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FIG. 3: Mass-Radius relation for magnetic white dwarfs for one Landau level systems, when (a) 
ej? = 2, (b) ep = 20. M is expressed in the units of Mq and R in the units of 10^ cm. The 
maximum mass is determined from the maximum value of pc for one level system as determined 
from equation Q. For ep = 2, the corresponding pc is 7.4 x 10^ gm/cc and for ep = 20, the 
corresponding pc is 1.2 x 10^^ gm/cc. 
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of ep. Thus the mass corresponding to that radius gives the maximum mass possible for the 
white dwarf. 

We know from the non-magnetic case that a stiffer equation of state can balance gravity 
more effectively leading to a higher value of radius (non-relativistic case). The same argu- 
ment can be applied here for F = 2.92 (lower density and hence outer region of the star), 
owing to the fact that the stiffness arises from the effect of magnetic field on degeneracy 
pressure. Hence from the mass-radius relation, mass would turn out to be higher than the 
usual one at a particular radius. Here we must be careful that even though the resulting 
analytical form of equation apparently does not have any restriction on the value of the 
maximum mass or radius, the maximum mass is restricted by the typical values of central 
density for white dwarf js] and the maximum Fermi energy of the star determined by equa- 
tion iQ. Figure |3] however shows that the magnetic white dwarfs have typically higher mass 
and higher radius. In general for typical densities of white dwarf, the mass of the white 
dwarf can exceed the Chandrasekhar mass limit. For ep = 2, however, we see from the 
mass-radius relation that the typical radius turns out to be higher than that of ep = 20, but 
the Chandrasekhar mass limit is not exceeded. 

In the equation of state for one level, F decreases from 2.92 in the low density region to 
the value 2 in the high density region. Quantitatively for a two level system, we can see 
from the equation of state that after the transition from the ground level to the first level, 
the slope of the curve changes drastically to very low values and then rises again. At this 
transition point, the pressure is nearly independent of density. This is physically improbable 
and the corresponding mass-radius relation would give an unstable branch. If a star starts 
out in this density range, it cannot form a compact star and will go into a runaway process. 



VI. DISCUSSION AND SUMMARY 



In this work we have calculated the equation of state of a degenerate electron gas in high 
magnetic field at zero temperature analytically. For simplicity, and keeping it analytically 
tractable, we have discussed the extreme cases when there are only one and two Landau levels 
for the system and studied the mass-radius relation for the case of one Landau level. Our 
pure analytical approach helps in understanding the underlying physics in great details. In 
the forthcoming work, we will show that the detailed numerical unrestricted solutions indeed 
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match with our analytical results 15 1. 

To make the problem theoretically tractable easily, we have chosen constant magnetic 
field throughout. However, this does not matter for the present purpose due to the following 
reason. Although the matter density and the magnetic field both vary from center (highest) 
to surface (lowest), by the time the density falls to about half the value of the central density, 
the mass generally crosses the Chandrasekhar mass limit. Thus, although we have considered 
a constant magnetic field, the field strength plays its major role and brings in new results in 
the high density regions of the star only. Hence it appears equivalent to the central magnetic 
field of the star. Now following previous work IJ, ll6[, one can adopt an inhomogeneous 
magnetic field profile in any compact star which is nearly constant throughout most of the 
star and then gradually falls off close to the surface (see Figure 5(b) in jl^). Thus the 
choice of an inhomogeneous magnetic field would not affect our main finding. A detailed 
description for the same will be given in a follow up paper based on detailed numerical 
analyses |l5[ |. 

Note that the size of a white dwarf is very large so that the effect of general relativity 
is insignificant there. In addition, the central density of the star is almost an order of 
magnitude larger than the energy density arised due to the corresponding magnetic field. 
As the effect of magnetic field arises in the high density regime only when the matter density 
exceeds the energy density arised due to the magnetic field, even the magnetic field does not 
contribute to the gravitational field leaving the system Newtonian. 

It is found that the ground level of the system is occupied first and then the higher levels 
start getting filled in. This transition appears as a kink in the pressure- density plot, what 
we have addressed analytically. For one level system that we have studied explicitly for the 
mass-radius relation (our numerical solutions addressed in a separate paper [155 would show 
the mass-radius relations for multi level systems), it is known to have no kink, for two level 
systems one kink, and so on for higher level systems, what we have studied analytically. 
We have derived a general expression for the positions of the kinks in the equation of state 
in terms of the magnetic field and other constants and the relation among them. It has 
been found that they bear a simple ratio, thus if we know exactly the value of density for 
a kink, then we can find out the density for other kinks. One of the cases considered here 
is for the Fermi energy 20meC^ and the system is of a one level. The minimum magnetic 
field required for this system is found to be 8.7 x 10^^ G. The corresponding mass of the 



17 



resulting white dwarf has been found to be larger than that predicted by Chandrasekhar 
and is about 2.3Mq when the radius is 8 x 10^ cm. We have also analysed analytically how 
is the mass-radius relation dependent on the Fermi energy and the magnetic field. 

We end by addressing the possible reason for not observing such a high field yet in a 
white dwarf. This could be due to the magnetic screening effects on the surface of the star. 
If the white dwarf is an accreting one, then the current in accreting plasma depositing on 
the surface of the white dwarf could create an induced magnetic moment of sign opposite to 
that of the original magnetic dipole, thus reducing the surface magnetic field of the white 
dwarf unaffecting the central field. In addition, the surface field could be several orders 
of magnitude smaller than the field in the central region which brings in the main results. 
Hence by estimating the surface field one should not interpret the rest. 
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